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Abstract 

We discuss the constraints on the inhomogeneous reheating scenario. In particular, 
we discuss the prospects for low scale inflation with a Hubble constant of the order 
of the gravitino mass, and the possibility that an MSSM flat direction is responsible 
for the density fluctuations. Thermal effects are generically fatal for the scenario, 
and can only be avoided for small enough mass and couplings of the flat direction 
field, excluding MSSM flat directions. Prompt decay at the end of inflation bypasses 
thermal constraints, and is compatible with both low scale inflation and MSSM flat 
directions. However, the level of non-Gaussianity is acceptably small only for a small 
window of inflaton mass and couplings. The upshot is that tuning of parameters is 
needed for the inhomogeneous reheating scenario to work. 



1 Introduction 



The inflationary paradigm in which the inflaton sector is responsible for the density 
perturbations is economical. One potential serves many purposes: driving inflation, 
ending inflation, and generating the observed spectrum of density perturbations. 
The resulting models of inflation are predictive, and falsifiable by experiment. The 
flip side of the coin, though, is that they are restrictive too, and it has proven 
extremely hard to build realistic models of inflation. Often a considerable amount 
of fine-tuning is needed for the model to satisfy all constraints. For this reason it 
may be worthwhile to explore alternatives. 

Inflationary models can be "liberated" if its task list is reduced |T]. This is 
the idea behind both the curvaton scenario [21 E] and the inhomogeneous reheating 
scenario Q , in which not the inflaton field but some other field is responsible for the 
density perturbations. The inflaton sector merely serves to drive and end inflation, 
and is considerably less constraint. Of course, the price to be paid is that a new 
field has to be introduced into the theory. But the costs are minimized if the new 
field and/or the new scales introduced are already present in our models of particle 
physics. There is already an extensive literature exploring this possibility in the 
context of the curvature scenario jSJ IH] , but little has been done in the context of the 
inhomogeneous reheating scenario (but see [3 El)- In this paper we will concentrate 
on the inhomogeneous reheating scenario. In particular, we will address the following 
two questions. Is low scale inflation — with the Hubble scale during inflation of the 
order of the gravitino mass — possible, such that the inflaton sector can be naturally 
identified with SUSY breaking sector? Can any of the flat directions in the minimal 
super symmetric standard model (MSSM) be responsible for the density fluctuations? 

Low scale inflation is hard to realize in the conventional setting in which the infla- 
ton is responsible for the density fluctuations. ^ In that case the density fluctuations 
are proportional to H^/Mpi, with the subscript * denoting the quantity evaluated 
at the time observable scales leave the horizon, some 60 e-folds before the end of 
inflation. Unless the slow roll parameter e = {Mp\V' /V)"^ <^ 1, which often requires 
some amount of fine tuning, high scale inflation is needed to get the observed size 
of density perturbations. In the curvaton scenario the density perturbations are 
proportional to H^/cx^, with a the vacuum expectation value (VEV) of the curvaton 
field. Since the VEV can be much lower than the Planck scale, low scale inflation is 
more natural in the curvaton scenario. However, in JU] it was shown that for the 
scenario to work the scale of inflation cannot be too low: > 10'' GeV. 

In the inhomogeneous reheating scenario the density perturbations are propor- 
tional to H/S^ or H/M, depending on the form of the decay rate, where S is the 
VEV of the flat direction field responsible for the density perturbations and M is 
some cutoff scale. Since both M and can be much smaller than the Planck scale, 

^The exception are "new inflation" models, in which the inflaton starts out close to the origin 
in]. When Taylor expanded around the inflaton VEV, the linear term is anomolous small, making 
this different behaviour possible. 
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also here, low scale inflation seems natural. We will see, that indeed, low scale infla- 
tion with of the order of the gravitino mass is possible in this context, although 
only in a small part of parameter space. The fluctuating decay rate scenario can 
then be married with any of the interesting inflaton models, in which the inflaton 
originates from the SUSY breaking sector [11]. 

In both the curvaton and inhomogeneous reheating scenario a scalar fleld other 
than the inflaton is responsible for the density perturbations. This scalar has to be 
light compared with the scale of inflation to be able to fluctuate freely, and produce 
the scale invariant perturbation spectrum observed. Obvious candidates for this 
scalar field are the MSSM flat directions. Within the curvaton scenario this possi- 
bility has been studied in [F. It was found that only under very special conditions 
can the MSSM scalar produce the density perturbations. The main problem is that 
the curvaton has to (nearly) dominate the energy density before decay. This means 
large initial VEVs and thus potential problems with non-renormalizable operators, 
and this means late decay and thus potential problems with disastrous finite tem- 
perature effects. In this respect, the inhomogeneous reheating scenario offers better 
prospects, as the MSSM flat direction does not need to dominate the energy density 
at the time of decay. Consequently, its VEV can be small enough to avoid prob- 
lems with the non-renormalizable potential, and the density perturbations may be 
generated before thermal effects become important. 

This paper is organized as follows. In the next section, we provide the back 
ground material. We describe those features of the inflaton and MSSM sector im- 
portant for the inhomogeneous reheating scenario. In section IHl we turn to the 
description of the various constraints on the inhomogeneous reheating scenario. Sec- 
tion 0] discusses the various time scales in the problem, giving further insights in the 
nature of especially the thermal constraints. We address the prospects for model 
building, concentrating on the possibility of low scale inflation and using MSSM flat 
directions, for a polynomial decay rate in section |3[ while section El discusses the 
same issues in the context of an (approximate) constant decay rate. In section [7| we 
discuss the varying mass scenario, which is a variation of the inhomogeneous reheat- 
ing scenario, and point out its virtues and drawbacks for model building. Finally, 
we conclude in section |H1 

2 Preliminaries 

The idea behind the inhomogeneous reheating scenario is the following. Consider 
the decay of the inflaton fleld, or more generally the fleld dominating the energy 
density (FDED), into standard model degrees of freedom. Suppose now that the 
decay rate for this process depends on the vacuum expectation value of some fleld. 
In supersymmetric theories as well as in superstring inspired theories, the effective 
couplings are functions of the various flelds in the theory, and thus so is the decay 
rate. If one of these flelds, call it S*, is light during inflation it can condense with 
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a large VEV. Moreover, it fluctuates freely during inflation. As a result, the decay 
rate is spatially fluctuating on superhorizon scales. 

FDED decay will happen at slightly different times in different parts of the 
universe. The regions in which decay has taken place are filled with radiation, 
while the not yet decayed regions are matter dominated. The universe expands at a 
different rate in different regions, resulting in fluctuations in the reheat temperature, 
hence, to adiabatic density perturbations. 

2.1 Decay rates 

The decay rate of the fleld dominating the energy density can be schematically 
written as F = X'^mK, with m the FDED mass, A its coupling, and K a phase 
space factor. All three quantities can have a field dependence, leading to different 
realizations of the inhomogeneous reheating scenario. The various possibilities will 
be discussed in detail in later sections. For now we just want to remark that we can 
divide the decay rates into two general classes, which differ in their iS-dependence. 
Here and in the following S denotes the ffat direction ( "flaton" ) fleld responsible for 
the density perturbations. 

The first class consists of decay rates polynomial in S, of the form 



Decay rates of this form can arise from non-normalizable operators. If there are 
several fields Si that can give rise to a decay rate of the above form, the dominant 
contribution will come from the field with the largest VEV. A large VEV favors 
directions with a small mass. This may naturally select a fiaton with ms < 0.1 if*, 
leading to a scale invariant spectrum. 

The second class of decay rates are of the form 



with S < M. This is the p — limit of the polynomial decay rate. At zeroth 
order the decay rate is constant, the S dependence comes only in through higher 
order terms. Such decay rates may be obtained from normalizable operators, or 
from phase space effects. We will refer to the above decay rates as respectively 
"polynomial" and (approximately) "constant". 

2.2 Inflaton sector 

We will not be concerned with the speciflc origin of inflation. All that is needed is 
that there is period of inflation, long enough to solve the horizon and flatness prob- 
lem, and with a Hubble constant that is almost constant \H^/H^\ -C 1. However, 
the mass of the fleld that stores most of the energy at the end of inflation will 
be an important parameter, so we digress somewhat on that. 



(p>i), 



(1) 





(2) 
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In one-field models of inflation the inflaton mass is bounded by the Hubble scale 
during inflation: < H^. Since the density perturbations are not produced by the 
inflaton field, the slow roll parameter 77 does not need to be much smaller than unity. 
Fast roll inflation is possible for > H^, but the number of e- folds is negligible 
small unless 

There is more freedom in multiple field models of inflation. Consider for example 
hybrid inflation with a potential 

^(0, x) = Vo + 5V{x) - \m\<^^ + \Xx^^ + (3) 

In supersymmetric theories the couplings A' and A are related. If the inflaton x is 
responsible for the density perturbations, A' has to be sufficiently small to assure a 
scale invariant spectrum, but this requirement can be relaxed in the inhomogeneous 
reheating scenario. Inflation occurs in the regime X > Xc = ^4>/V\'i and the 
inflaton rolls slowly in the potential 5V{x)- Inflation ends when x < Xci and 
acquires a VEV 

where we have set = in the vacuum. The field dominating the energy density 
at the end of inflation is the waterfall field 0. Taking A ~ 1, the mass of this field 
in the true vacuum is ~ V^^'^ and thus ^ ~ Vq. ^ Another example of 
two-field inflation is the recently proposed "new old inflation" J^I- The potential is 
of the hybrid inflation type, and also in this case ^ H^. is possible. 

In the next sections the mass denotes the mass of the field dominating the 
energy density at the end of inflation, i.e., the inflaton field in one field models of 
inflation and the waterfall field in hybrid inflation. With an abuse of language will 
refer to this field in both cases as the inflaton field. We parametrize 

~ (IFI), 

~/5v/H:, (2FI), (5) 

with /3 < 1 for both one field inflation (IFI) and two field inflation (2FI) such as 
hybrid inflation. Since it is not expected that ^ Vq, the /3-factor is maximum 
for 2FI. 

A variation of the inhomogeneous reheating scenario was proposed in pj. Sup- 
pose that the inflaton decays into heavy particles if), whose mass is set by a flat 
direction VEV. The particles freeze-out, and when they become non-relativistic 
they soon come to dominate the energy density in the universe. The density per- 
turbations are generated during the decay of the ip particles, and are sourced in this 

^We will use units in which the reduced Planck mass Afpi = SttG — 1 
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case by a varying mass as opposed to a varying coupling. ^ This scenario seems 
less economical since it needs the introduction of yet another field. We consider it 
for completeness though. The mass is bounded by the temperature < T 
for thermal production. If they are direct infiaton decay products, < for 
perturbative decay, and < 10^ — lO^m,^ for bosons and fermions respectively in 
non-perturbative decay 

We will only consider perturbative infiaton decay. It is well known that the 
waterfall field in hybrid infiation generically decays non-perturbatively, in a rapid 
process dubbed instant preheating [15]. The negative (mass)^ at the end of infiation 
leads to a spinodal instability, and the condensate breaks up almost instantly. All 
long wavelength modes below some critical wavelength k^, are excited. The waterfall 
field still dominates the energy density, but now the energy is not only stored in 
the zero-mode but in all modes with k < k^. If decay happens before annihilation 
reactions become important, the inhomogeneous reheating scenario still works as 
before. 

2.3 Low scale inflation 

Several infiationary models have been constructed in which the infiaton sector is 
linked to SUSY breaking ^1]. Then naturally Vq^^ ~ 10~* and ~ Vq^'^ ~ "^3/2 ~ 
10~^^ for gravity mediated SUSY breaking. In anomaly mediated SUSY breaking 
schemes, the scale of infiation Vq^^ is larger by one or two orders of magnitude, 
whereas for gauge mediation the scale is smaller. 

2.4 MSSM flat directions 

The fiat directions of the MSSM consist of gauge invariant operators composed 
of MSSM scalar fields p!B|. This polynomial is commonly parametrized as X = 
0", with n the dimension of X. A non-zero VEV for (p will break the standard 
model gauge symmetry. All fields entering in the fiat directions which are left after 
the Higgs mechanism (except the linear combination which receives the VEV after 
diagonalization) have masses of order h(j), with h a gauge coupling, due to their D- 
term couplings to the fiat direction VEV jT7j. In addition superpotential couplings 
of the form W = h(j)qq lead to effective mass terms Mq ~ h(j) with h the MSSM 
Yukawas. 

Therefore, if the fiaton 5* responsible for the density perturbations is identified 
with one of the MSSM fiatons 0, at least one of its couplings is of gauge strength 
h 0.1. The Yukawas vary between 10~^ < h < 1. 

■^Of course, one could also consider a decay rate for ^p with a varying coupling constant. This 
is the same as considering infiaton decay with a varying coupling, only less economical, since an 
extra field is introduced. 
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Figure 1: Efficiency parameter a as a function of T^/Hend 

3 Constraints 

We list here the various constraints on the inhomogeneous reheating scenario. 
3.1 Density perturbations 

The density perturbations can be parametrized by the gauge invariant quantity (, 
which describes the the density perturbations on uniform curvature shces JHl- The 
perturbation generated by a fluctuating inflaton decay rate are [U [71 [TU] 

C = c^, (6) 

with a the efficiency parameter. is to be evaluated at the time of decay, when 
H ~ r^. After reheating is completed the metric can be written as ds"^ = —dt^ + 
g'^{V^)tdx^ j2nj. The function g oc F^" parametrizes the difference in expansion rates 
in different regions, and thus the resulting density perturbations. It can be calculated 
numerically by integrating the coupled equations of motion for the inflaton and the 
radiation bath, together with the Friedman equation. The slope gives the efficiency 
parameter a, which is shown in Fig. 1. 

The parameter a depends on the ratio F^/iJcnd with Hend the Hubble constant 
at the end of inflation. If F^ ^ -ffend, decay occurs almost instantaneously in all 
regions of space, almost independently of the fluctuations in the decay rate, and the 
resulting density perturbations will be small. In the opposite limit F^ -C -ffend the 
efficiency is maximal, and a reaches its maximum value a = 1/6, which agrees with 
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analytical estimates 01[7inni- We can parametrize 



1/6, < lO-^i/end, 

{ 1/6 - 1/15, < 10-2 - lifend, (7) 

0.1(ifend/r</,), > i^end- 

The value for ( measured by WMAP is C = 2 x 10"^ [21 . The right amount of 
density perturbations for a polynomial decay rate of the form Eq. are obtained 
if 



6S_ 



_ 2 X 10 

dec ap 



-5 



(8) 



The equivalent for the constant decay rate of Eq. Q is 



S'^-^SS 



Ml 



(9) 



dec aqr 



where we have used S'/M -C 1. 

If the density perturbations are produced during the decay by a field ip other 
than the infiaton field, that dominates the energy density sometime after the end of 
inflation, then generically -C i^end and a = 1/6. If the FDED has a varying mass, 
density perturbations are already produced before decay, as domination happens at 
different times in different parts of the universe. The resulting perturbations after 
decay are 

6 1^ S 

to be evaluated at the time of decay. If the mass oc S, the decay rate is 
polynomial and the density perturbations are up to a minus sign given by Eq. (jH)) 
with p = 1 and a = 1/6. 

In the standard scenario in which the infiaton is responsible for the density 
perturbations the non-Gaussianity is small. The slow roll conditions assure that 
the perturbations are Gaussian at production, and as the adiabatic perturbations 
remain constant on superhorizon scales, no further non-Gaussianity is produced. In 
contrast, large non-Gaussianity is possible in the inhomogeneous reheating scenario. 
Apart from the possibility of non-Gaussianity at production, non-Gaussianities can 
be produced due to the non-linear evolution if S starts rolling in the potential, due 
to the non-linear relation between the decay rate and 6S, and due to the non-linear 
relation between the decay rate and the final perturbations |71 1201 [22] • The general 
rule is that the less efficient the transfer of 5* perturbations in metric perturbations 
is, the larger the non-Gaussianities. This is easy to understand, since if the transfer 
is less efficient, larger fluctuations SS/S are needed to obtain the observed spectrum. 



^While finishing up this paper \22\ appeared, claiming the density perturbations produced in 
the fluctuating mass case is about a factor 10 larger than the estimate Eq. 10), which is taken 
from 7_ . The difference will not affect our conclusions in an essential way. 
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A level of Gaussianity as probed by WMAP is assured if 



SS < lOS, fill 



both at the time of production as well as at decay. In particularly, this means 
55^, ~ H^:/{2'k) < 105**. The transfer of perturbations is inefficient if the transfer 
parameter a is small, or for the constant decay rate if the ratio S/M is small. 
Indeed, combining Eqs. (jHUHl) with the constraint from Gaussianity, Eq. (jllj) . bounds 
a > 10"^ and a{Sdec/My > 10"^ respectively. When V""(^*) < the flaton zero 
mode and its evolution starts rolling in its potential. This may introduce additional 
inefficiencies, which will be discussed in the next subsection. 

The produced spectrum of perturbations is nearly scale invariant, as dictated by 
observations, if the mass of the fluctuating field is sufficiently small ms ^ O.IH^. 
Note that the same must hold for the effective mass mes = V" generated by higher 
order terms in the potential. A quartic term V = k\S\^ will lift the flatness of the 
potential, unless k is exceedingly small. For this reason we will not consider such 
potentials in the following. One should also consider non-renormalizable operators 
in the potential of the form 

_NW^ (,2) 

with A some cutoff scale, e.g. the Planck mass. In an expanding universe supersym- 
metry is broken dynamically, leading to soft mass terms for all scalars of the order 
of the Hubble constant. Such large masses m ~ if will spoil the scale invariance of 
the perturbations, and should be avoided. We will assume that the Hubble induced 
flaton mass is at least one order of magnitude below its canonical value, assuring it 
is unimportant at all times. Soft mass terms can be suppressed by either invoking 
symmetries, or by allowing tuning. 

The density perturbations produced trough the inhomogeneous reheating mech- 
anism are by assumption the dominant ones; other sources should lead to negligi- 
ble perturbations. In particular, the perturbations produced by the inflaton fleld 
should be negligible small, which is assured for if* < iimax = 10^^. Further, at 
the time of S'-decay ps/ptotai ^ 1, otherwise the flaton perturbations will dominate 
through the curvaton mechanism. In the curvaton scenario ( ~ r/{27r)H^/ S^, with 
r = ps/ptotai ~ Slg~^ and g the largest flaton coupling, determining its decay rate. 
The curvaton perturbations are sub-dominant if g^S^. 



3.2 Evolution of S and 6S 

The decay rate is a function of the flat direction fleld, and therefore a function of 
time. The inflaton decays when H/T^{S) drops below unity. This excludes cases in 
which the decay rate decreases faster than the Hubble constant. 

Once the effective flaton mass becomes of the order of the Hubble rate, V"{S^) ~ 
H^, the flaton starts oscillating in its potential. For a quadratic potential ps ~ a~^. 
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During inflaton oscillations, the scale factor of the universe scales as a oc H^'^^^, 
and S (X H. The ratio H/r^{S) remains constant, and the inflaton does not decay. 
If the potential is dominated by non-renormalizable operators, then S rolls down 
approximately critically damped V^'j^ ~ H^. We parametrize the post-inflationary 
evolution of the flaton field 

S ^ fS,, (13) 

with 

/ oc < _±_ (14) 

where = (l/2)m||S'p and Vnr given by Eq. (IT^ . 

Further, the inflaton should decay before the decay of the flat direction field: 
> Ts- This is automatically satisfied if > ms and flaton couplings less 

than unity. The same constraints hold if it is not the inflaton but some other field 

dominating the energy density, which has a varying decay rate. 

As the Hubble constant drops below ^yV"{S) not only the zero mode but also 

the fluctuations 5S start evolving. The equations of motion for the zero mode and 

the superhorizon fluctuations are 

S + 3HS + V'{S) = 0, 

6S + 3H6S + V"{S)6S = 0. (15) 

For a quadratic potential, the equations of motion for S and 6S are identical (to 
linear order) and the ratio {6S/S) remains constant. If the potential is dominated 
by non-renormalizable operators the fluctuations will be damped with respect to the 
zero-mode according to fIE\ 

n— 4 



D OC i/2{n-2)^ (16) 

with the damping factor D defined as 



s J \s , ^ 

There is no damping for n = 4 operators. 

Using the definitions of / and D we can express the equations for the density 
perturbations Eqs. (jHl El) in terms of the quantities during inflation: 

S^^Sx 10^ apDH,, (18) 

for a polynomial decay rate. Here we have used 55** ~ H^/{27t). Gaussianity as in 
Eq. (fTT|) requires D > 10~^p~^((l/6)/a). The equivalent expression for a constant 
decay rate is 

^ « (19) 

Ml aqrDfi ^ ' 

with D/" > 10-3(^M/^,)9((l/6)/a)(gr)-i to assure Gaussianity. 

Apart from the zero temperature potential the evolution of S and 5S can be 
sourced by finite temperature effects. 
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3.3 Thermal constraints 



The radiation bath affects the ffaton condensate in two ways, through thermal scat- 
tering and through thermal corrections to the ffaton mass. The fields that couple to 
the ffat directions through a gauge/Yukawa coupling h, denote them by x, have a 
large effective mass ~ hS. The x particles are in thermal equilibrium with the 
radiation bath, if their effective mass is smaller than the temperature 

m^^hS< T. (20) 

It is important to realize that even before inffaton decay has completed, there is a 
dilute plasma with temperature j2Sl 

T = [TlHf/^ = {V^HY'\ (21) 

with Tr ^ a/T^ the reheat temperature at the end of the reheating process when 
r<^ ~ H. In between the end of inffation and inffaton decay the effective temper- 
ature scales as T oc H^/^^ instead of the T oc H^/'^ during radiation domination. 
The plasma reaches its maximum temperature Tmax immediately after the end of 
inffation: ^ 

T^ax ~ ~ (i^end^<^)^/^ (22) 

with i^end the Hubble constant at the end of inffation. 

The back reaction of the x quanta on the ffat direction field induces thermal 
corrections to the ffaton mass. If the x quanta are in thermal equilibrium, when 
Eq. ()20|) is satisfied, there is an induced "plasma mass" of the form 

Sm% ~ /i^^ ~ h^T\ {hS < T), (23) 

where the last equality hold if the x quanta have a thermal distribution with ra^ ~ 
and ~ T. There is also a thermal correction in the opposite limit — T < m^, 
and the x quanta are out of equilibrium — as a result of integrating out the heavy 
degrees of freedom |2Z||2H]- The massless vector and chiral superfields, i.e., those not 
"Higgsed" by the ffat direction VEV S, generate a two-loop free energy proportional 
to h''^T'^, with h' the corresponding gauge or Yukawa coupling. The running of h' 
changes at the renormalization scale /i ~ ~ hS, if is charged under the 
gauge group, respectively couples to the matter field in question. Integrating out 
the heavy x fields above this scale, leads to an effective "RG potential" of the form: 

SVnc CTh"THog (^^^ , ihS>T), (24) 

with ct a constant of order one which can have either positive or negative sign, 
depending on the matter content. The sign is negative if the running is dominated 
by integrating out Higgsed gauge bosons. 
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Large thermal masses, rrif^ > H > ms, induce early oscillations of the flat 
direction field. During these oscillation, the energy density stored in the fiat direction 
scales as 1211 

Ps{H) = ^^,%Ps{Ho). (25) 
m{Ho) 

Thermal evaporation of the fiat direction through collisions with the x particles 
occurs if (1) the scattering rate is large Fgcat ^ H, (2) the x particles are in thermal 
equilibrium, and (3) the energy density stored in the thermal bath is larger than the 
energy density in the fiat direction. This last condition will generically be satisfied. 
The equihbrium condition is given by Eq. ((201) • The scattering rate is Fscat ~ ^x'^- 
For a superpotential term W = hSxx fhe cross section for xS scattering is a ~ 
h'^a/E^^, where we have assumed scattering is dominated by fermions. The typical 
center of mass energy is Ecm ~ y/Tm, the mean of the typical x energy (~ T) 
and S energy (~ m = hT, where the last equality holds when thermal masses are 
important). Then Fscat ~ hg^T, with ~ 0.1 a gauge coupling [77j . 

Early thermal evaporation obviously kills the inhomogeneous reheating scenario. 
As we will see, early induced oscillations are generically also fatal. Therefore, for 
the fluctuating decay scenario to work, thermal effects should be negligible. There 
are three possibilities. The first is that the particles coupling to the fiat directions 
are out of equilibrium. Thermal scattering is negligible. There is a renormalization 
group (RG) induced potential of the form Eq (j24|) . If the thermal mass is less than 
the Hubble scale SV^q <^ H^, the fields remains frozen and thermal effects are 
negligible. 

The second possibility is that the effective plasmon mass is lower than the tem- 
perature. Thermal scattering is delayed with respect to the onset of induced oscil- 
lations. Therefore, thermal effects are absent it the induced plasma mass, Eq (|23|). 
is smaller than the Hubble rate. 

The time when thermal effects become important may be delayed if the thermal 
plasma is initially far from a thermal equilibrium distribution. If n^/E^ -C T^, the 
plasma mass is then below its equilibrium value. The scattering rate may also be 
suppressed, as it scales with n^. ^ Likewise, if the x s-i'e non-relativistic, the RG 
potential can be suppressed if the free energy of the massless field are below their 
equilibrium value. 

An initial distribution with n^^/E^ <^ can happen if the plasmons are direct 
infiaton decay products with E^ ~ ^ T. Or else, if the plasmons are not in- 
flaton decay products, their initial number density is small. The plasmons acquire 
chemical equilibrium when the rate for number changing interactions Fint becomes 
of the order of the Hubble constant. Fint ~ h'^i9^T, where hpi is the coupling of x 
to the thermal bath and g a gauge coupling. If the plasmons carry gauge charges, 
^pi ~ fl' ~ 10~^; the time scale for chemical equilibrium Fint ~ 10~'^T is smaller than 

^Note that there is an enhancement in the scattering rate for smaller induced mass, as cr cx 

I/Wth- 
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the plasma mass /iT, and the onset of plasma induced oscillations can be delayed 
only for h > 10~^. 

The induced thermal masses and scattering rates are a function of the temper- 
ature T ~ {HT fj,{S)Y^^ . Hence, the constraints depend on how the inflaton decay 
rate depends on S, and are therefore model dependent. The dependence of the 
plasma temperature on if, Eq. (j2ip . is derived for a constant inflaton decay rate. 
We will approximate the temperature by this formula also in the case of a variable 
decay rate, so that the temperature becomes a function of both H and S. 

If oscillations of the fiat direction field are set off by thermally induced masses, 
the field amplitude decreases as given by Eq. (j^Hl)- For plasma masses ~ hT 
this implies 

/0,^7/(8+p)^ (26) 

with p = for a constant decay rate. Since the potential is approximately quadratic 
D ~ 1. We postpone a discussion of / and D for RG induced potentials to the next 
section. 

To avoid gravitino overproduction and thereby spoiling the successful nucleosyn- 
thesis predictions, the reheat temperature has to be sufficiently small Tr < 10~^°, or 
equivalently V < 10~^^ [30j. This is not a hard bound, in the sense that depending 
on the specifics of the SUSY breaking mechanism, the bound can shift some orders 
of magnitude. Moreover, the gravitino problem can be solved for example by invok- 
ing a period of thermal inflation [3T] . A bound that cannot be tampered with is the 
big bang nucleosynthesis (BBN) bound: The field dominating the energy density 
should decay before BBN, Tr > 10^22 or T > lO^^^ 

In this context, we should also remark that the baryon asymmetry of the universe 
has to be created either during or after FDED decay. If FDED decay happens below 
the electroweak scale < 10~^^, sphalerons are out of equilibrium, and an Affieck- 
Dine like mechanism for baryogenesis is required. 



4 Time scales 

The important time scales in the problem are the following. Inflation ends at 

^^end ~ C^*, (27) 

with ( parameterizing the difference between the time observable scales leave the 
horizon, and the end of inflation. For simplicity we set C ~ 1, although lower values 
are possible. For example, in IFI with a quartic potential ( ~ 0.1. 

The large VEV of the inflaton gives a large mass to its decay products, and 
inflaton decay is kinematically forbidden during inflation. Decay becomes possible 
as soon as the inflaton starts oscillating in its potential. Therefore, the Hubble 
constant at the time of decay is H^ec ~ niin[r(^, i^end]- 
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The scales at which the vacuum mass and effective mass generated by non 
renormahzable operators become important are respectively H ms and 



NR 



(28) 



We denote the scale at which thermal effects become important by iJth- The 
particles coupling to the ffat direction are in equilibrium for T > hS, or equivalently 



H > i^eq, with 



~ min 



end 



(29) 



Then for H > ifgq, -f^th ~ Hpi the time scale at which the plasma mass of Eq. ()23|) 



becomes of the order of the Hubble scale, 



while for H < H^q, Hth 



scale at which the RG potential of Eq. ()24j) becomes important. Here 



pi 



RG 



mm 



mm 



end 
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2p + 1 + p{p — l)log 



hS . 

rj~\ J 



end 



i^RG the time 

(30) 
(31) 



where p = corresponds to the constant decay rate. As discussed in section El 
evaporation is delayed with respect to the onset of thermally induced oscillations. 
The onset of thermal effects can be delayed if the initial plasma is far from an 
equilibrium distribution, and the number densities of the particles giving mass to 
S is less than the equilibrium value <^ T^. The delay is until number changing 
interactions become important, at ifdei ~ Tint with 



Hdci ~ min 



(v/)^/^^y^ if, 



end 



(32) 



with hp\ the coupling of the particle responsible for the induced thermal mass of S, 
i.e., the plasmon x which couples to S if the x's are relativistic, and the light degrees 
of freedom whose free energy is altered by the running of x for non-relativistic x's- 



We can then distinguish the following cases: 

1. > -fiend- Infiaton decay is kinematically inaccessible during inflation, and 
the infiaton decays promptly at the end of inflation. 

2. The infiaton decays after the end of inflation, but while the flaton is still frozen, 
i.e., before the onset of oscillations of the flat direction field. Consequently, 
there is no damping D = 1 and S'dec = S*. This is the case for 

-f^end > > ms, Hth, -f^NR- (33) 

3. Infiaton decay occurs after the end of inflation, and after the fiaton field has 
started rolling in the quadratic potential, the non-renormalizable potential 
and/or the thermally induced potential. Damping and the evolution of S 
should be taken into account, and / and D are generically different from 
unity. 
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Case 1 Case 1 — prompt decay — is obtained if > i/cnd- The transfer of 
flaton perturbations to metric perturbations is inefficient in this hmit, and large 
non-Gaussianity is possible. Gaussianity as in Eq. (fTT|) and the requirement of 
prompt decay together constrain the decay rate, or equivalently the parameter Tq 
(see Eqs. dUE))), to the range 

{5pxlO'H,)-P< ^ <5pxW'H;P, (r^pol.) (34) 

1< ^ ^ 5gr X 10^ (^^y , (r.cnst.) (35) 

for a polynomial and constant decay rate respectively. Here we have used Eqs. (ll8[fT9|) 
There is only a small window of Tq values for which prompt decay is possible. 

There are no thermal or damping effects, and the only constraints come from 
the observed magnitude of perturbations Eqs. (fTHl IT^ with f = D = 1, and scale 
invariance m^, a/^r(5'*) < O.IH^. In addition the gravitino problem has to be 
addressed for H^, > 10~^^. 



Case 2 Inflaton decay occurs after the end of inflation, but before the Hubble 
constant drops below the effective mass of the fiat direction field, which aside 
from the vacuum mass can have contributions from thermal effects and from non- 
renormalizable operators in the Lagrangian. The zero-mode and fluctuations are 
frozen, and f = D = 1. 

The flaton vacuum mass ms, respectively the effective mass from non-renormalizable 
operators, is smaller than the Hubble constant at inflaton decay for 

Hdec > ms, -ffNR- (36) 

Thermal effects play no role if the plasmons are in equilibrium, but inflaton decay 
happens before the plasma mass becomes important, that is if 

Hdec > Hpi, Heq (37) 

It should be understood that Hpi,Heq < -ffend, i-e., plasma masses are not impor- 
tant immediately after inflation has ended, nor are the plasmons out of equilibrium 
from the start. Another possibility for the thermal effects to be absent is that the 
plasmons fall out of equilibrium before the equilibrium thermal effects become im- 
portant. This includes the case in which the plasmons are out of equilibrium right 
from the start. A RG correction to the potential is generated, which should be 
negligible H^ig < -^dec- Both conditions are satisfled for 

h < S''\ (38) 

Finally, thermal effects may not set in until H ~ i/dei, if the plasma is initially far 
from the thermal equilibrium distribution. Thermal effects are avoided if 

-f^del ^ Hdec (39) 
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Case 3 The flaton starts rolling in its potential before inflaton decay. A polyno- 
mial decay rate changes with changing S, whereas a constant decay rate remains 
(approximately) constant. Hence, these two cases are quite different and we will 
discuss them separately. 

As said, the decay rate is time-dependent for a polynomial decay rate. The first 
remark to be made is that the expression for the density perturbations Eq. (jH)) is 
derived for a constant decay rate. However, we expect that the time-dependence 
will not significantly alter the size of the final perturbations. 

Since the decay rate changes with time, there is the possibility that the decay 
rate decreases faster than the Hubble constant, and inflaton decay will not occur. 
In fact, this is the case if the fiaton starts oscillating in the quadratic potential, 
when H < m^. The fiaton red shifts S oc H, see Eq. (jHj), and decay is impossible 
for all p > 1. Similarly, decay does not occur when the fiaton starts rolling in the 
non-renormalizable potential, and n < p + 2. If oscillations of the flat direction field 
are set off by a plasma mass rripi ~ hT, the field amplitude decreases as given by 
Eq. (j26j) . For p > 2 the decay rate decreases faster than the Hubble rate, and the 
inflaton cannot decay during thermally induced oscillations either. Therefore 

p>n+2 (40) 

should be satisfied for the inhomogeneous reheating scenario to be possible. Here 
the superscript eq indicates that the plasmons coupling to the flat direction are in 
thermal equilibrium, and the subscripts p > 2 and p > n + 2 indicate that only 
decay rates satisfying these conditions are implied. 

For the same Hubble constant during inflation, the flat direction VEV at inflaton 
decay is smaller (larger) than in case 1-2, due to both the damping of the fluctuations 
and the evolution of 5*. Now S = /S** = 8 x lO^apfDH^, with a ^ 1/6. For 
fD < 1 {FD > 1), inflaton decay occurs for Hubble constants which are a factor 
{fDy smaller (larger) compared to case 2. The bounds on the fiaton mass and 
couphngs, Eqs. IHTl IHHl Unj) . are stronger (weaker) by factors of {fDy. In 
addition, there is the new constraint that the damping should not be too large, 
D > 10~^, to assure a Gaussian perturbation spectrum. 

Consider first the effect of the non-renormalizable potential V^r of Eq. (fT^ with 
n > p + 2, so that inflaton decay is possible. If H^r > T^, the flat direction 
field starts slow rolling in the potential before the epoch of inflaton decay. For 
simplicity, we will concentrate on the case that thermal effects are negligible at all 
times -ffNR > > ms,Hth_. Then fD < 1 and all constraints are stronger. The 
decay rate at if < H^r is 

/ TT \2/(n-2) 

r^H) = U— r<^(^NR) (41) 

V-"NR/ 

Decay occurs at r0(ifdec) ~ -f^dec- The damping factor is D = (ifdecZ-f^NR)^"""^^^^^""^^- 
The amplitude 5* can be found as a function of by solving the equation 5** ~ 
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lO'^pDH^. The amplitude at the time of inflaton decay is Sdec ~ (-f^dec/ -^^nr)^^^" ^''•S'*. 
For example, for p = 2 and taking k = 1 and A = 1, this leads to 



(42) 



The RG induced potential becomes important before inflaton decay if 

i/eq > i^pl, i^dec, & /^>5l/'. (43) 

The behavior of the zero mode and the fluctuations depends critically on the sign 
of ct- The equations of motions are 

S + 3HS + AHS"-^ [p\og{^) + 1^ 5 = 0, 

5S + "iHS'S + AHS'P-^ {2P-I+ p{p - 1) log(^)^ 6S = 0, (44) 

with A = CT^h'^To, which can have either sign depending on the sign of cy. For 
positive sign ct > and p > 2 the effective potential is up to factors of log(5') given 
by Vth oc 5*^. The field rolls towards lower values, and / and D are approximately as 
given in Eqs. (fT^ ll(j |l . The thermal potential shuts off when hS < H, at which point 
a plasma mass hT is generated. If Hpi < H the field freezes, and as the temperature 
drops the plasmons again fall out of equilibrium. The induced RG potential leads to 
the decrease of S until the plasmons regain equilibrium. And so forth. As a result, 
the VEV tracks hS ~ T. Inflaton decay should happen before the flaton mass 
or plasma effects become important, i.e., Eqs. (jSHl ISZ|) should be satisfied. If the 
tracking is halted due to the effects of the non-renormalizable potential, S decreases 
further, but the inhomogeneous reheating scenario is still possible for n > p + 2. 
Both the effects of Vnr and 6V with ct > lead to fD < 1 and all bounds are 
stronger compared to case 2. 

The situation is completely different for c-r < 0. Now the potential is minimized 
for large S. When damping is negligible, for p >2 the instability leads to exponen- 
tial growth of the zero mode. For p > 3 the growth is stopped by damping, when 
i^RG ~ H^, and the VEV tracks S*^"^ ~ H/A. For p = 2 however, the thermal 
mass is up to log factors independent of S and damping remains unimportant. The 
growth of the zero mode is halted instead when non-renormalizable terms become 
important, or when the inflaton decays, whatever comes first. If the log term in 
Eq. ()44j) does not already dominate initially, with the the fast growth of S it will 
soon do so. Then for p = 2 the zero mode and the fluctuations roll in the same 
effective potential and their ratio remains constant, ~ 1. Since f'D ^ 1, inflaton 
decay is earlier making constraints weaker. Moreover, since S increases, the out of 
equilibrium condition hS > T remains valid. 
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For an approximately constant decay rate of the form Eq. Q the time evolution 
of the flaton will not alter the time of decay. Inflaton decay can still occur, for 
example, after the flaton starts oscillating in its quadratic potential. 

For the same Hubble constant during inflation, the flat direction VEV at inflaton 
decay is smaller than in case 1-2, due to both the damping of the fluctuations and the 
evolution of S, see Eq. (fT^ . The decay rate does not alter, but the bounds on H^, S** 
and M to get the right density perturbation are stronger for f^D < 1. Specifically, 
taking Gaussianity and the requirement S < M into account, constrains > 
10~^. For oscillations set off by the mass term, the non-renormalizable potential, or 
by plasma masses f^D < 1. From the expressions for / and D, Eqs. (fT^ 1161 OUj) . it 
follows that 

mc,, i/pi<10% & i7NR<10% (45) 

for the varying decay rate scenario to work. The flaton VEV and should be 
tuned even more than for cases 1 and 2. 

Things may be different for the RG potential of Eq. ()24|). For a constant decay 
rate, the zero mode and fluctuations have mass terms with opposite sign 

AH 

S + 3HS+--S = 0, 
AH 

S'S + 3HSS - —-6S = (46) 

with A = CT^h'^TQ. For ct > the zero mode decreases, while the fluctuations 
increase, and vise versa for c-r < 0. For c-r > 0, the decrease of S is halted when the 
plasmons fall out of equilibrium, and hS tracks T. The interesting case is g = 1, 
since only then f^D > 1. For ct > the growth of the zero-mode is halted by 
damping, and 6V" ~ H^. Only for g > 3 is fW oc i7-('2-2)/2 > i. 

However, there is not much that can be done to improve the situation. This can 
easily seen by looking at the expression for the density perturbations at the time of 
decay, Eq. Q, which implies SS^ec > 10~"^M. Gaussianity and the cutoff constrain 
lO^^dec < 'S'dec < M. For a given cutoff and Hubble scale during inflation, there is 
only a small window for S'dec- For ct > 0, S increases towards earlier time, but it 
cannot excess the cutoff scale. For ct < 0, S decreases towards earlier time whereas 
SS increases, leading to large non-Gaussianities when S ~ 6S. There is no room 
for S to either decrease or increase much, and a considerable amount of tuning is 
needed if the RG potential is to dominate at some point. 



5 Inflaton decay through non-renormalizable op- 
erators 

A polynomial decay rate, as in Eq. ([Q), can arise from non-renormalizable superpo- 
tential couplings of the form 4j 

W = Xo^HQk (47) 
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with M the cutoff scale, and S = (qi) 7^ is the field responsible for the density 
perturbations. This corresponds to 5-dim non-renormalizable operators in the po- 
tential. Such operators may be obtained from integrating out physics above the 
cutoff scale M. For infiaton decay into MSSM degrees of freedom the combination 
QiQjflk should form a gauge invariant. For S part of a MSSM flat direction, the only 
combinations are qqq, qqq and hqq with q and q left-handed quark/lepton superfields 
and their charge conjugate. Infiaton decay is unsuppressed for > nig, with nig 
the lightest quark/lepton fields in the superfields qj, q^. For S a SM singlet, the only 
possibility is Sh^hd ^ 

The effective coupling is A = \o{S/M). Then the decay rate is of the form Eq. (|T)) 
with ^ 

with m^i, = (3H^ {I3\/1TI) for IFI (2FI), and < 1. The value Fq is bounded from 
above by Aq, /? ^ 1. The lower bound comes from the requirement that the infiaton 
decays before BBN: 

10-47 

5.1 Low scale inflation 

Is the fluctuating decay scenario with a decay rate of the form Eq. (fTlHSj) compatible 
with low scale inflation ~ ""^3/2 ~ 10^^^? And if so, what are the bounds on the 
couplings and mass of the flaton field? 

The decay rate depends on the specifics of the infiaton sector. However, Fq can 
be bounded by BBN and the scale of inflation to lie in the range 10~^^/(/Da)^ < 
Fo ~ 10-^^(10-^°)/?A^/M2 for IFI (2FI). For M ~ 1 there is no parameter space for 
IFI. The cutoff should be larger than the flat direction VEV M > IQ-^'^aDH^. 
The gravitino problem is absent for Fq < 10^/(/D)^. 

Case 1 Prompt decay is possible for 10^° ^ To ^ 10^^, as follows from Eqs. (j34jl . 
Such large decay rates are not possible in IFI, while for 2FI it can be obtained for 
10^0 < (3Xl/M^ < 1029. This means 10"^^ < M < 10^^° and Ag/5 > 10"^ The 
scale M cannot be identified with the scale of inflation, and a new scale has to be 
introduced in the theory. 

The flat direction couplings are unconstrained, whereas its mass ms ^ 10~^^, in 
order to obtain a scale invariant spectrum. The effective mass generated by the non- 
renormalizable potential is also sufficiently small for Planck suppressed operators 
with K, A ~ 1. Note however, that if these operators are suppressed by the same 

^In this case a normalizable superpotential coupling W = \'(t>qjqk is also possible. The decay 
through non-renormalizable operators dominates for XqS/AI > A'. 

''Higher values of p can be obtained through interactions W = Xo{S/M)"^(j)qjqk with m > 2. 
Since the constraints are stronger for higher values of p we will not discuss this possibility. 
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scale as the effective inffation coupling, i.e., A ~ M, then n = 4 operators should 
be absent or suppressed k <^ 1. Further, M > 10~^^ (10~^^) for n = 5 (n = 6) 
operators. 

The gravitino problem has to be addressed. 

Case 2 For smaller decay rates, Fq < lO^^D'^, the inflaton decays well after the 
end of inflation and case 2 and 3 apply. There are additional constraints, of which 
the strongest ones comes from finite temperature effects: early thermally induced 
oscillations and thermal evaporation. The inflaton decays while the flat direction 
field is still frozen, and case 2 applies, if these thermal effects are negligible small. 

A word of caution here. All estimates are order of magnitude estimates. But 
this is even more so when thermal effects are considered. The reason is that the 
thermally induced mass by particles in (out of) equilibrium is a good approximation 
in the limit hS <^ T {hS ^ T). However, in the limit hS T both approximations 
break down. 

The first possibility for thermal effects to be absent is that the particles coupling 
to S fall out of equilibrium before plasma effects become important, and with the 
RG potential small Eq. (jHHjl . This is only possible for small rates Fq, inconsistent 
with the BBN bound. 

The second possibility is that the plasmons are in equilibrium until inflaton 
decay, and the plasma mass is small, Eq. (j37j) . This requires h < 10~^^Fy^. Taking 
M ~ 1, this implies very small couplings and mass: h < 10~^^ and ms ^ 10"'^^ for 
2FI. Identifying the cutoff with the scale of inflation M ~ \/H^ ~ 10~^ improves 
the situation considerably, but still h < IQ-^^ (10"^^) and ms < 10~20 (10"^^) for 
2FI (IFI). The constraints are weakest in the limit M 10^^^, but at the expense 
of introducing a new scale in the system, and of gravitino over production. 

Finally, there is the possibility that the initial plasma is far from an equilibrium 
distribution, and the onset of thermal effects is delayed, Eq. flHUj) . H^ei ^ F^ gives the 
same constraints as in the previous paragraph, with the replacement h — * lO'^/ipj. 
Here hpi is the coupling of the x particles (those particles coupling to S) to the 
plasma. If x has gauge charges hpi ~ 0.1, and delay of thermal effects is not 
possible. 

The non-renormalizable operators are sub-dominant until inflaton decay, and 
Eq. dnni) is satisfied, for Fq > 10-^^^''-^\k / A)""-^ . For M, A ~ 1 and 2FI n = 4 
operators should be suppressed. For M ~ 10^^ and A ~ 1 non-renormalizable 
operators are consistent with 2FI, and with IFI for n > 5. However, if M ~ A ~ 
10"^ one needs n > 5 (n > 6) for 2FI (IFI). 

Case 3 Case 3 applies if the flat direction field starts rolling in its potential before 
inflaton decay. If the potential is dominated by the mass term ms, the plasma mass 
mpi ~ hT, or n = 4 non-renormalizable operators the decay rate decreases faster 
than the Hubble constant, and decay cannot occur. Non-renormalizable operators 
with n > 5, and RG induced potentials with ct > have fD < 1 due to evolution of 
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the zero-mode and its fluctuations. The aheady tight constraints of case 2 become 
stronger by appropriate factors of fD. 

The only possibly interesting case is when the potential is dominated by RG 
effects with ct < 0, see Eq. ()24|) . Both the zero- mode and its fluctuations grow 
exponentially, so that / ^ 1 and D ~ 1, and hence fD ^1. A RG potential is 
generated when (see Eq. (jl^ ) 

h > min[10i3/2 V?;, 10'/\ron (50) 

and therefore Fq < 10~^° for couplings less than one. Moreover, ii h > S^^^ ~ 10-^3/2 
the induced thermal mass is larger than the decay width, and the inflaton start 
rolling in SVth before decay. The small rates Tq required are naturally for M ~ 1. 
A large cutoff has the additional advantage that the zero mode can grow by a huge 
amount without exceeding the cutoff. Note also, that small Tq is consistent with 
the BBN bound for / > 1. 

For all Fq values under consideration ^ , the zero-mode starts growing rapidly as 
soon as if ~ -f^Rc ~ To/?-^, until its growth is halted by either the non-renormalizable 
potential (5Vrg ~ Vnr), or inflaton decay (F^ ~ H). Approximating the in- 
crease to be instantaneous, the maximum value of the zero mode before the non- 
renormalizable potential becomes important is ~ -f^RG ~ Fq/i^. Decay happens 
before this maximum value is reached if 

r^>i^RG =^ ^max>/i' (51) 

with = FoS'^ax- Eq. (jSOI) implies large couplings, whereas Eq. (j^H) is satisfied 
more easily for small couplings. This contradiction can only be remedied for large 
n and small decay rates. Decay before Vnr becomes important does not happen 
for n = 4, 5 operators; it is marginally consistent with n = 6 operators in the 
limit Fo ^ 10-26 and h ^ IQ-^^/^. The flaton mass has to be extremely small: 
ms < F^ ~ Hkg ~ 10-3^. 

One could contemplate the possibility that the growth of 5* is halted by non- 
renormalizable operators, and the flaton slow rolls in the non-renormalizable poten- 
tial before decay. But this is to no avail. In this case the left hand side of Eq. 
is F^ ~ (/-D)NR'S'max' with (/-D)nr oc H parameterizing the evolution of S and 6S 
in the non-renormalizable potential. The right hand side is likewise proportional 
to H. Thus for H < -f^NR? both sides decrease with the same rate, and it remains 
impossible to satisfy Eq. (j^D). In other words, if decay does not happen before the 
non-renormalizable potential becomes important, it will certainly not happen dur- 
ing the slow roll in the non-renormalizable potential if the right amount of density 
perturbations are to be produced. 

^For 10"^^ < To < 10^^" the growth of the zero mode can start immediately at the end of 
inflation, at H ^ Hcnd- The non-renormahzable potential becomes important before decay for all 
n. 
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5.2 MSSM flat directions 



Can the flat directions of tlie MSSM play the role of S for a decay rate of the form 
Eq. PHj) ? We take ms ~ "^.3/2 ~ 10~^^. The Yukawa couplings of the MSSM vary 
between 10~^ and 1, whereas the gauge couplings h ~ 0.1. As discussed in the 
previous section, only if the inflaton decays promptly at the end of inflation, is low 
scale inflation consistent with MSSM flatons. But what are the conditions for an 
arbitrary scale of inflation in the range lOm^ < < 10~^? 

Since the inflaton decay rate has to be greater than the flat direction mass there 
is always a gravitino problem. 

Case 1 Prompt decay occurs for 10^® < TqH^ < 10^, see Eq. fl34|) . This is not 
possible for IFI for any scale of inflation. For 2FI prompt decay is possible for 
> 10ms, M < lO^H*^^ < 1, and 5* < M, which is hardest to satisfy in the limit 

^ i^max ~ 10-^ 

Low scale inflation, with the Hubble constant of the order of the gravitino mass, 
comes at the cost of introducing a new scale M which cannot be identified with 
any of the known scales, such as the SUSY breaking scale, the GUT scale, or the 
Planck scale. The scale M can be identified with the scale of inflation, M ~ \/H^ 
for > 10-'/i(3'Xl). 

Case 2 For smaller Fq the inflaton decays well after the end of inflation, and 
cases 2 and 3 apply. The decay rate is greater than the MSSM soft mass for 
H* ^ 10^^'^(/-D)^^rQ Requiring < 10^^ and decay after the end of infla- 
tion constrains lO~^H^/{fDY < Fq < 10^^/(D^if*), where we have taken a ^ 1/6. 

Let's first consider case 2. There are three ways to avoid thermally induced 
early oscillations. The first possibility is that the plasmons coupling to S fall out 
of equilibrium before plasma effects become important, and the non-equilibrium 
thermal mass is small Eq. (j38|) . However, for MSSM gauge couplings the non- 
equilibrium thermal mass is never small, and this possibility is excluded. 

The second possibility for the thermal effects to be negligible is that the plasmons 
are in equilibrium until inflaton decay, and the plasma mass is small, Eq. ()37p . 
Plasma masses are negligible for h < y/TT^, which excludes MSSM gauge couplings. 

The third possibility is that the thermal plasma is initially far from an equilib- 
rium distribution. However, for gauge particles the thermal effects can be delayed 
at most until H^ei ^ O.IH^. This offers not much perspective either. 

Case 3 Just as for low scale inflation, all thermal constraints are stronger for 
case 3 if fD < 1. Hence, case 3 is likewise incompatible with MSSM flatons. The 
only possible exception are MSSM flat directions whose potential is dominated by a 
tachyonic RG mass, i.e., Eq. ()24|1 with ct < 0, so that fD > 1 is possible. The zero- 
mode grows exponentially once the tachyonic mass exceeds the Hubble constant. As 
discussed in the previous subsection, the inhomogeneous reheating scenario can only 
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work if the inflaton decays before the growth is halted by the non-renormahzable 
potential, and Eq. is satisfied. 

Oscillations induced by the plasma mass reduce the flaton amplitude. As a result 
Heq oc H^^^^, and the plasmons remain in equilibrium until decay. The only way 
out is that (some of the) the particles coupling to 5* are out of equilibrium from the 
start ifeq > ^^end, or 

To < (52) 

The MSSM flaton couples to several fields. If one of these couplings h' is small so 
that the above equation is not satisfied, then in addition i^Rc > -f^* and i^RG > Hpi 
for the RG induced mass to dominate over the plasma mass from the start. This 
is only possible for flat direction with top Yukawa interactions /i ~ 1; in all other 
cases Eq (jHSl) should be satisfied for all flaton couplings. 

Consider first MSSM fiat directions with top Yukawas. The RG induced potential 
can dominate over all other contributions to the potential immediately after the end 
of inflation. The zero mode and fluctuations grow exponentially. Decay happens 
before the non-renormalizable potential becomes important for H^, < Tq < lO^iJ* 
and To > 10~^, 10~^^ for n = 5, 6 operators; it is not possible for n = 4 operators. 

For flat directions without top Yukawa couplings the growth of the zero-mode 
starts at H^q ~ h^T^ < H^. Moreover, Eq. should be satisfied. These con- 
straints together are severe, prohibiting decay before V^-r becomes important for 
n < 5, whereas it is only marginally allowed for n = 6 in the limit Fq ~ i/o ^ 10~^. 



6 Constant decay rate 

Renormalizable couplings Consider a superpotential coupling of the form ^ 

W = MHuHd (53) 



with the coupling 



A = An 



+ 



(54) 



with the ellipses denoting higher order terms. Since the higher order corrections play 
no role we will omit them in the following. In supersymmetric and string inspired 
models, all couplings and masses are, rather than being constants, functions of 
scalar fields in the theory. The higher order corrections can also arise from non- 
renormalizable operators, such as those in Eq. ()47|) . Within the MSSM the above 
coupling is the only gauge invariant possibility. 
The decay rate is 



1 + 1^ 

' M 



H 2 



To = 



(55) 
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which is of the form Eq. Q with g > 1, r = 2 and M some cutoff scale. As before 
= (3H^ for IFI (2FI) with /9 < 1. The decay rate is bounded from below 

by BBN 

To > 10-^^ (56) 

In addition the mass of the inflaton should be larger than that of the decay products, 
otherwise decay is kinematically forbidden. 

The density perturbations are given by Eq. (fTU|) . All scales in the problem have 
to lie close together, < 10-^M/(a/«D) and < < M. The amount of fine- 
tuning is increased by inefficiencies, when af^D < 1/6. Considering the density 
perturbations in terms of the variable at the time of decay, Eq. 0, similarly leads 
to the conclusion that the amount of fine tuning is increased also in the opposite 
limit pD > 1. 

Phase space factor Consider 2-body decay of the infiaton, through a coupling 
of the form 

W = \o(t)^^ (57) 

The decay rate is 




2m^y _ 



To = ^m^, (58) 

where the square root comes from integration over phase space. We will refer to 
this decay rate as the phase space (PS) decay rate. If the mass of ip is set through 
a coupling to a fiat direction, i.e., = XS, the decay rate is of the form Eq. 
with g = 2, r = l/2 and M = m^/{2h). Now M is not a fundamental scale in 
the problem. The density perturbations are given by Eq. (fT^ . Also in this case 
M ~ 10 - lO^H^. This gives the relation 

/ lO-lO^A, (IFI), 
lO-lO^Av^, (2FI). 

The lower bound on Fq comes from BBN, Eq. (I56j) . There is also an upper 
bound, from Aq < 1 and /3 < 1 (lO^y/K) for IFI (2FI). 

Note that the back reaction of ip on the fiat direction, induces a thermal mass 
for S. Said in another way, the couplings of the fiaton field to which we refer as h, 
and which control the strength of the thermal effects, include also A. 



6.1 Low scale inflation 

Can the inhomogeneous reheating scenario work for low scale infiation H^, ~ ^3/2 ~ 
10^^^ for decay rates of the form Eqs. ((21 ESI EH))? To obtain the right density 
fiuctuations requires S** ~ 1 — lO^if* and M ~ 10 — 10'^ if* for all examples considered 
above. The less efficient the mechanism — a{S^/MY small — the closer the scales 
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lie together. M cannot be identified with inflationary scale, and thus introduces a 
new scale in the problem. This is not the case for the PS decay rate, where M is 
related to the effective mass of the inflaton decay products, and does not represent 
a fundamental scale. 

The decay rate is bounded by the BBN constraint and the scale of inflation 
10"^^ < To ~ 10-^8(10-i°)A§/3 for IFI (2FI). For the PS decay rate (3 < (3^^^ ~ IQ-^ 
for 2FI as follows from Eq. (j59|) . There is a gravitino problem for Fq > 10~^^. 

Case 1 Prompt decay is possible for 10~^^ < Fq < 10^^^{H^/My. Small values 
for q are favored. Note that for the upper limit, which saturates the Gaussianity 
constraint, ~ and thus H^/M < 0.1. Prompt decay is not possible for 
IFI, whereas for 2FI it constrains lO'^ < p\l < IQ^^iHjMf. It is also is also 
incompatible with a PS decay rate. 

The window of allowed inflaton mass and couplings is much smaller than for 
the polynomial decay rate. One reason is that the perturbations are transferred 
less efficiently to the radiation bath, due to the {S^/MY suppression factors. The 
second reason is that the constraint 1 < ^<f>/H < 10^ confines /3Ao also within 
three decades. This is in contrast with the polynomial decay rate, where due to the 
S, and thus a, dependence of F<^, the parameter combination /SAq/M^ is constraint 
only within nine decades. 

The other flat direction couplings are unconstrained, whereas its mass < 
10~^^ for a scale invariant spectrum. The effective mass from the non-renormalizable 
potential is also sufficiently small for Planck suppressed operators with k, A ~ 1. 
However, if these operators are generated at the same scale as the effective cutoff, 
A ~ M, then the non-renormalizable operators are only marginally consistent for 
all n. 

The gravitino problem has to be addressed. 

Case 2 For smaller decay rates, Fq < 10^^^, decay occurs well after the end of 
inflation, which brings us to case 2 and 3. Non-renormalizable operators play no role 
for Fo > 10-28k/A, IO-^^k/A^, IQ-^^k/A^ for n = 4, 5, 6. If k, A ~ 1, then only 
n = 4 operators need to be considered. Non-renormalizable operators with A ~ M 
always dominate before decay if Fq < 10~^^. 

Let's consider case 2 flrst. The thermal constraints are severe. They can be sat- 
isfled if the particles coupling to S are out of equilibrium, while the non-equilibrium 
thermal potential remains unimportant, see Eq. (jHHj) . The non-equilibrium potential 
is only small for couplings X, h < 10~^. But such small couplings are inconsistent 
with the plasmons being out of equilibrium. 

The second way in which thermal effects can be avoided, is that the plasmons are 
in equilibrium, and the plasma mass is small, see Eq. ()37|) . The plasma mass is only 
negligible for couplings X, h < vTq. The couplings scale with v^Fq, and must be 
especially small for IFF The flaton mass is bounded by ms < Fq- The constraints 
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are weakest in the limit Tq — > 10^^^, i.e., when decay happens shortly after the end 
of inflation. 

Finally, the thermal effects are delayed if the initial plasma is far from an equi- 
librium distribution. The ensuing constraints are the same as in the previous para- 
graph with the replacement h —>■ 10~'^h'^Y- The plasma couplings have to be small 

/ipi < io-3(ro/io-i6)i/4. 

Case 3 The non-renormalizable potential generically plays no role. If either the 
zero temperature mass or the plasma mass becomes important before infiaton decay, 
this implies fine-tuning the already tuned values of iJ*, 5*, M. And without opening 
parameter space much, see Eq. (jl^ . 

The RG potential offers no better prospects, whether ct is positive or negative. 
The VEV S is bounded by the scale of inflation and the cutoff M, both at the 
end of inflation and at the time of infiaton decay. To obtain the observed density 
perturbations, all scales have to lie within four decades of each other at all times, 
and consequently there is not much room for evolution of 5* and 6S. 

6.2 MSSM flat direction 

Can the MSSM flat directions be responsible for the density fluctuations, for a 
constant decay rate of the form Eqs. (j21 ESI MSSM scalars have soft SUSY 

breaking mass ms ~ 10~^^, and at least one or more gauge couplings with /i ~ 0.1. 
What are the conditions on and Tq for a successful scenario? 

Case 1 Just as for low scale inflation, prompt decay is only possible for a 2FI 
within a small window lO'^y/lT^ < (3X1 < 10^(if*/M)*. The scale of inflation can 
be identified with the cutoff scale only for ~ 10~^ and M ~ \^H7 ~ 10^^. In 
this case the infiaton mass and coupling have to be large i3Xq — > 1. Prompt decay 
is incompatible with a PS decay rate for all scales of inflation. 

The flaton coupling and mass are unconstrained, except for ms < O.liJ* from 
scale invariance. The effective mass generated by non-renormalizable operators is 
likewise sufficiently small for cutoffs A > M. 

The gravitino problem has to be addressed. 

Case 2 There are two possibilities for thermal effects to be absent. Either the 
induced plasma mass is negligible, or if the plasmons are out of equilibrium, the 
induced RG mass is sub-dominant, Eqs. IHHj) . This is the case for h < \/T^, 
respectively h < ~ 10^/H^. Since if* < 10^^, and Fq < (otherwise prompt 
decay), both options are inconsistent with gauge couplings of order h ~ 0.1. 

The onset of thermal effects can be delayed, see Eq. (jHOI), only for 10~^, 
and even then F<^ > O.lii*. Hence, this hardly opens up parameter space. 
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Case 3 If decay happens well after the end of inflation, thermal effects are always 
important. Decay should follow quickly after the onset of thermally induced motion 
of S and 6S. This is because Gaussianity and the cutoff constrain SS, S, and M 
to lie close to each other both at the end of inflation and at the time of decay, and 
there is little room for evolution of S and 6S. Therefore, ifth ^ lO^F,^ or equivalently 
Fq ^ 10~^Hj,. Only a small bit of parameter space is opened up, and only at the 
cost of fine tuning all parameters. 



7 Varying mass scenario 

Suppose the field dominating the energy density has an effective mass set by a 
couphng to a flat direction ~ AS*, with A a gauge/Yukawa coupling. If the ip 
field couples to SM degrees of freedom with a coupling Aq the decay width is 

r, ~ ^xs. (60) 

Decay is unsuppressed if is larger than the masses of the decay products. The 
decay width is polynomial in S, of the form Eq. ([[]), with 

ro = ^, P=i. (61) 

Stt 

The observed perturbations are obtained, if Eq. (fTSj) is satisfied with p = 1. The 
decay rate decreases with S, and ms < E^ for decay to take place at all. The 
thermal bath originates from inflaton decay, and thus the thermal time scales are 
as given in Eqs. IHUl EH) with E^ the ^-independent inflaton decay rate. 

There are important differences with the examples discussed up till now, in which 
the inflaton was the EDED. Now, a new field is introduced in the theory, that is 
to dominate the energy density sometime after inflaton decay. The new field is 
accompanied by new parameters, making the model less predictive. There is no 
equivalent of case 1, as the inhomogeneous reheating mechanism only works if ip 
decays after it (nearly) dominates the energy density, E^ < -ffdom- The thermal 
bath is produced by inflaton decay, and in the period between -ffend < H < T^, 
thermal effects should be taken into account. Thus, in contrast with the inflaton as 
EDED, a thermal bath is already present before ip domination. 



7.1 Constraints 

Here we discuss the conditions specific to the varying mass scenario, namely the 
requirement of domination, and the thermal effects before domination. 

Consider first the case in which the ip quanta are initially in thermal equilibrium. 
Eurther assume that the annihilation rate is smaller than Eann ^ "o^^- ^ Then 

^For larger annihilation rates freeze-out occurs at lower temperature T ^ /20, and the 
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the ijj particles come to dominate the energy density as soon as they become non- 
relativistic, at Tdom ~ or -ffdom ~ As was first noted in [7j, the ip plasmons 
induce a thermal mass for the fiat direction field 5* of the form ms ~ AT. Requiring 
this thermal mass to be sub-dominant at all times gives 

mth < i^dcc < i^dom ^ S>Mpi (62) 

Such large VEVs are hard to reconcile with low scale infiation if* < i^max, and with 
the presence of non-renormalizable operators. What is more, for an initial trans- 
Planckian VEV the fiat direction field itself will come to dominate the energy density 
of the universe while still frozen, leading to an ^-dominated period of infiation. 

The assumptions made in the above argument is that the ■^/'-quanta reach thermal 
equilibrium before domination, and that the FDED decays while S is still frozen. 
If either one of these assumptions is dropped, the negative conclusion implied by 
Eq. (jU^ may be avoided. We discuss some possibilities in turn. 

A) One possibility is that the ip quanta are direct infiaton decay products, which 
are out of equilibrium from the beginning 

i^i>^T ^/^^ (63) 

This requires large infiaton masses > m^. The initial distribution is non-thermal, 
and n^/s will remain constant if in addition the annihilation rate small. This means 
sufficiently small ip couplings, ruling out the possibility of identifying ip with an 
MSSM field. Moreover the couplings between ip and the MSSM sector have to be 
sufficiently small so that tjj decays after domination: 

Tv, < ^dom- (64) 

An RG potential of the form Eq. ()24j) is induced, which is negligible only for H^ig < 
or 

'"-^{jj ■ ^^^^ 

Once again small couplings are needed, ruling out the possibility of identifying S 
with a MSSM fiat direction. Otherwise, if the RG potential does get important 
and ct > 0, this will lead to damping, leading to non-Gaussianity for D < 10~^. 
For Ct < on the other hand, S decreases until the ■j/'-quanta acquire thermal 
equilibrium. This brings us straight back to the constraint Eq. (I62j) . 

To summarize, small couplings are needed for the annihilation rate and the 
induced thermal mass to be small, and MSSM scalars cannot play a role. There is 
a tension with the out of equilibrium condition, which requires a large coupling. 

B) The ip field dominates the energy density immediately after the end of infia- 
tion. This requires the infiaton to decay (almost) exclusively into ip quanta. Note 

number density is Boltzmann suppressed. Domination then happens at Hubble constants 
Hd ^ and the constraint in Eq. (|62|l becomes much stronger 
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that although Eq. (jU^ is trivially avoided in this way, soon after the end of inflation 
either the plasma mass or the RG potential generated by the coupling W = \Sipip 
will become important. The situation is the same as for a polynomial decay rate, 
as discussed in section V, except that case 1 is excluded. 

C) The inflaton decays non-perturbatively, through resonance effects. The initial 
distribution of ip particles is far from equilibrium, as well as all the other particles in 
the universe. The highly non-linear, non-equilibrium character of preheating makes 
it impossible to make any definite predictions. 

D) The rate of number changing interactions is small, and chemical equilibrium 
is not attained. The initial distribution is non-thermal if the ^/''s are direct inflaton 
decay products. The plasma mass m^i ~ hn^/E^ can be smaller than its equilibrium 
value. However, to avoid that -ffdom oc n^E^ is likewise smaller, and Eq. is not 
ameliorated, it should be E^ exceeding its equilibrium value rather than being 
below its equilibrium value. 

E) The thermal mass mpi ~ AT does become important before inflaton decay. 
The flaton field starts oscillating in the potential well with decreasing amplitude 
S oc if''/^, see Eq. (j2Sl)- If decays during the thermally induced oscillations, i.e., 
before the zero temperature potential becomes important, then / ~ (r^(S'^,)/ifpi)''/^ 
and D ~ 1. Unless r^(5'*) Hpi decay is delayed, and the thermal constraints, as 
well as the constraint < F^, are much stronger. 

7.2 Model building 

The varying mass scenario, in which the FDED is not the inflaton but some other 
field, is more elaborate than scenarios in which the FDED is the inflaton. There is 
an extra step: the production and subsequent domination of quanta. Not only 
does the introduction of an extra field lead to more parameters, and thus to less 
predictability, it also introduces extra constraints. In particular, there is already 
a thermal plasma before the domination of ip, and thermal effects should be taken 
into account. 

If the initial distribution of x quanta is an equilibrium one, the scenario does not 
work, because of Eq. (jHH). The natural way out is to assume the ip quanta never are 
in thermal equilibrium. However, non-equilibrium thermal effects still play a role, 
and it requires small coupling and/or tuning to make it work. 

The fluctuating mass scenario is not advantageous for low scale inflation. Prompt 
decay at the end of inflation is not possible, as decay should occur only after 
the ^/'-quanta come to dominate the universe. Late decay implies strong thermal 
constraints, and only very small flaton masses and couplings are consistent. The 
prospects are much worse than for the inflaton as FDED. 

Thermal constraints, which already play a role before domination, can generically 
only be avoided for sufficiently small couplings. This excludes the possibility of 
identifying ip with an MSSM scalar, as well as the possibility of identifying the 
flaton S with an MSSM flat direction. 
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8 Conclusions 



In the inhomogeneous reheating scenario not the inflaton but some other field is 
responsible for the observed density perturbations. It would be economical if the 
new fields and scales introduced in this scenario could be identified with the fields and 
scales appearing in our models of particle physics. In particular, we have addressed 
the following two questions. Is low scale inflation possible, with the Hubble scale of 
the order of the gravitino, such that the inflaton sector can be naturally identifled 
with the SUSY breaking sector? Can any of the MSSM flat directions be responsible 
for the density fluctuations? 

We discussed various decay rates, obtained from non-renormalizable couplings, 
renormalizable couplings, and from phase space effects. For the last two examples, 
the observed density perturbations can only be obtained if the Hubble constant, the 
flaton VEV, and the cutoff scale all lie within 4 decades. This requires some tuning, 
especially since Gaussianity of the perturbations is only assured for > lOH^. For 
the phase space decay rate the cutoff is not a fundamental scale in the theory. This 
has the advantage that there is no need to explain the origin of this scale. The 
disadvantage is that this model is quite constraint, and for example, prompt decay 
is excluded. For the non-renormalizable decay rate there is more freedom, since the 
density perturbations do not depend on the cutoff. 

After inflation, already before inflaton decay has completed, there is a dilute 
plasma. Fields coupling to the flat direction fleld, whether they are in equilibrium 
with the thermal bath or not, will induce a thermal mass for S. Such thermal cor- 
rections will lead to early induced oscillations of the flat direction fleld, which are 
generically fatal for the inhomogeneous reheating scenario. The thermal effects can 
be avoided if all flat direction couplings are small, which excludes MSSM flat di- 
rections. For a polynomial decay rate, which can originate from non-renormalizable 
couplings, in addition the flaton mass has to be much smaller than the Hubble scale 
during inflation. 

The thermal constraints can be trivially avoided if the inflaton decays promptly 
at the end of inflation, so that there is no time for thermal effects to act. In this case 
the flat direction mass and couplings are unconstrained, except for the requirement 
that nis ^ O.IH^. The reheat temperature is high, and there is a potential gravitino 
problem. Prompt decay however, is only consistent with a Gaussian perturbation 
spectrum for a small window of inflaton masses and couplings. The constraints 
are stronger for a constant decay rate. In particular, decay rates which inherit 
their S dependence from phase space effects are incompatible with prompt decay. 
Large inflaton mass and couplings are needed, inconsistent with one fleld models of 
inflation. 

Finally, we discussed a variation of the inhomogeneous reheating scenario, in 
which not the inflaton field but another field has a fiuctuating decay rate. After 
inflation, this field first has to come to dominate the energy density, and then decay. 
The prospects for model building are bleak. There is no analog of prompt decay: 
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thermal effects are always there, even before the field comes to dominate the energy 
density. Only with small couplings can one avoid the disastrous consequences of 
thermally induced masses. This excludes the possibility of identifying either the 
decaying field or the fiat direction field with an MSSM field. Moreover, there is an 
extra step in the model, making it less advantageous for low scale inflation. 
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